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We study idealizations of the full nonlinear Schwinger-Dyson equations for 
the asymptotically free theory of </> 3 in six dimensions in its meta-stable 
vacuum. We begin with the cubic nonlinearity and go on to all-order nonlin- 
earities which contain instanton effects. In an asymptotically free theory the 
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and perturbative methods fail from the outset. We show how our toy models 



43 of the cubic Schwinger-Dyson equations contain the usual diseases of pertur- 

Q_i' bation theory in the massless limit (e.g., factorially-divergent (3- functions, 



singular Borel-transform kernels associated with infrared renormalons) and 
show how these models yield specific mechanisms for removing such singulari- 
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ties when there is a mass gap. The solutions to these homogeneous equations, 
in spite of being ultraviolet finite, still depend on an undetermined param- 
eter equivalent to the perturbative renormalization scale /u. In the all-order 
nonlinear equation we show how to recover the usual renormalization-group- 
improved instanton effects and associated factorial divergences. 
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I. INTRODUCTION 



Asymptotically free theories appear to have perturbatively calculable properties in the ultravi- 
olet (UV) where the running coupling gets small. The well-known price paid for this convenience 
of perturbation theory in the UV is its complete breakdown at some low momentum scale, usually 
thought of as being of O(Arg) where Arg is the renormalization group (RG) mass. However, as 
one goes to higher orders A" in perturbation theory, the critical momentum scale grows exponen- 
tially in A" essentially because the number of graphs grows as N\. That is to say, the contribution 
of the A^ th order term in perturbation theory behaves like N\(ag 2 ) N where g ~ (In k 2 )~ l l 2 is the 
running charge; to keep this term small as A" grows requires an exponential increase in k 2 . There 
are also factorial divergence associated with renormalons JT[|. Consequently, perturbation theory 
to all orders cannot be used for any momentum scale, however large, in an asymptotically free 
theory without understanding how to deal with non-Borel-summable factorial divergences. This 
is currently not a practical limit on perturbation theory in QCD, where the critical momentum 
does not creep into the UV until A" is rather large, say AT > 10. 

A more practical issue is understanding QCD processes at infrared (IR) scales k ~ Arg and 
the attendant phenomena of confinement, condensates, renormalons, large instantons, etc. [0 The 
purpose of the present paper is to discuss the factorial divergences mentioned above, as well as IR 
issues (except for confinement) in a toy model of asymptotic freedom [^]: 3 in six dimensions (0g). 
This theory is terminally ill because the Hamiltonian is unbounded below but we can go quite far 
before encountering pathologies. The idea is to study the theory at all momentum scales from the 
IR to the UV in a connected way, without using the conventional crutches of perturbation theory. 
Our study will use the machinery of nonlinear Schwinger- Dyson equations j|] for 

Even without the spin complications of QCD, it is too hard for us to solve the full Schwinger- 
Dyson equations numerically, even if only the cubic nonlinearity is saved in the Schwinger-Dyson 
equation for the vertex. So we turn to toy models, first saving only the cubic term for the 
vertex and later considering some aspects of the all-order nonlinearities in the vertex equation. 
Although we have no proof, we believe that these toy models fairly represent qualitatively, but 
not quantitatively, all the diseases and their cures we would encounter in the full Schwinger-Dyson 
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equations. 

Of course, any study of 0g must ultimately break down because the theory itself is ill-defined; it 
has no stable vacuum. We will see that the breakdown occurs when we try to define the sum of the 
non-Borel-summable N\ divergences associated with vertex skeleton graphs having 2N + 1 vertices. 
These sums contain the contributions of instantons which we need not add as an explicit 

ingredient in analyzing the field theory; they are already in the Schwinger-Dyson equations. 

The Schwinger-Dyson equations of an asymptotically free theory have a special feature: the 
renormalization constants which appear in the renormalized Schwinger-Dyson equations must 
vanish. This vanishing is clear from the canonical form of the 0g vertex equation (or its analog in 
QCD), in which the vertex renormalization constant Z\ appears only as an additive inhomogeneous 
term. At large momentum, asymptotically freedom requires the vertex and all the other terms 
of the Schwinger-Dyson equation to vanish, as we will show; this is inconsistent unless Z\ = 0. 
However, one must be careful. The unqualified statement Z\ = is false and leads to paradoxes. 
A more precise statement resembles that of lattice gauge theory where the theory is first defined 
with an UV cutoff Auv and Z\ — > at a specific rate as Auv — > oo. 

We will see that it is correct simply to set Z\ = in the canonical form of the vertex equation 
(or in an analogous equation for the running coupling, which is RG-invariant). The resulting 
vertex equation is homogeneous, which completely disconnects it from any perturbative approach 
and vastly complicates the analysis in some ways. 

At this point we turn to toy models of the homogeneous and nonlinear vertex Schwinger-Dyson 
equation. Not all models that we have studied will be presented, since they all show the same 
features. We believe that all these features will occur in the full Schwinger-Dyson equation as 
well. Our studies (both analytic and numeric) fall into several categories: 

1 A- Cubic nonlinearity: massless. The toy-model vertex Schwinger-Dyson equation for the 
running coupling can be converted to a (nonlinear) ordinary differential equation. There are 
several types of solutions, all of them having singularities at sufficiently small momenta and having 
undetermined parameters. One of them behaves as expected from perturbation theory for large 
momentum where the running coupling varies as (In A; 2 ) -1 / 2 ; other solutions vanish as an inverse 
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power of k 2 (modulo logarithms). We find a differential equation for the /3-function whose power- 
series solution diverges like N\, with all terms negative 0; this divergence is associated with IR 
renormalons. We use the Laplace transform of the vertex equation in In k 2 to find an equation 
for the Borel transform of, e.g., the /5-function, and find explicitly its pole structure at the IR 
renormalon singularity. 

IB-Cubic nonlinearity: massive. In this case, as is well-known, all the IR renormalon difficul- 
ties are resolved (at least for large enough mass), and the toy-model Schwinger- Dyson equations 
are non-singular. The only solutions which are regular for all Euclidean momenta behave like 
(lnfc 2 ) -1 ' 2 for large k. All solutions depend on a single continuous real parameter, equivalent to 
the RG renormalization point /i; this is somewhat of a surprise, since all models yield UV-finite 
solutions needing no cutoff or other regularization. Modulo this parameter, the solutions show 
uniquely-determined power-law corrections to perturbation theory (e.g., condensates and other 
higher- twist terms). Generally, the toy models (as well as the full Schwinger-Dyson equation) can- 
not be reduced to differential equations, but there is one notable exception, which we analyze in 
some detail. As for the massless case, this is a cubically- nonlinear second-order ordinary differen- 
tial equation. The Laplace-transform form of the vertex equation shows that the renormalon pole 
singularity is cancelled by another term vanishing like an inverse power of k 2 (modulo logarithms) 
at large k. 

2. All-order nonlinearity. First, we prove a theorem (similar to another used recently |7| to 
study N\ divergences in 0|) that the imaginary part, in Minkowski space, of every skeleton graph 
for the 0g vertex has neither IR nor UV logarithms. The implication is that the UV and IR 
behavior of the dressed graphs is entirely determined by the dressed vertex itself, which allows 
for an approach to the all-order problem by successive approximations, beginning by inputting 
the solutions to the cubic Schwinger-Dyson equation. When we use the (ln/c 2 ) _1//2 behavior found 
earlier as input, we find that the iV-vertex graph has the asymptotic UV behavior ~ (In k 2 )^^ N ^ 2 ^ 2 . 
Using some deep graph-theoretic results and other tools developed in an earlier (f>\ study, we show 
that this momentum factor is multiplied by a N N\, where the positive constant a is rather close 
to what would be expected from a one-loop RG-improved instanton analysis. The sum over all 
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iV is ambiguous, but should have an imaginary part reflecting the instability of <f>\. The real part 
shows the usual power-law behavior expected from instantons. 

We comment on some of these results in more detail. Among our major results are our demon- 
strations of i) how the Schwinger- Dyson equation can be expressed as a nonlinear integral equa- 
tion for what would be called the Borel transform in perturbation theory, ii) how the massless 
Schwinger-Dyson equation gives rise to poles in this transform and iii) how the massive Schwinger- 
Dyson equation gets rid of these poles by condensate terms. In its simplest terms, the massless 
case yields a Borel integral for the /3-function derived from the running coupling of the form 

/3(g) = constant x f° da (1) 
Jo 2 — a 

where H(a) is regular at a = 2. This corresponds to a /3-function with terms behaving like 
— g 2N+1 (b/2) N N\ at large N. Here b is the lowest-order coefficient in the /3-function: (3 = —bg 3 +- ■ ■. 
We show that, when masses are included, Eq. [l| changes to 

(3(g) = constant x H da ( e -/(^ 2 ) _ e ~V^)\ (2) 

Jo 2 — a v ' 

which has no singularity at a = 2. The running coupling has a similar expression, but with 1 / (bg 2 ) 

replaced by ln(/c 2 + M 2 ); this shows that the cancellation involves terms vanishing like (k 2 )~ 2 at 

large k. Of course, the perturbative expansion of the /5-function is the same as before, but Eq. § 

shows how that sum is actually defined by the Schwinger-Dyson equation. 

Another interesting result is that all our toy models, while perfectly finite in the UV and 

requiring no cutoffs or regularizations, still have a single real parameter in the solutions, which is 

not fixed by the vertex Schwinger-Dyson equation above. This parameter is equivalent to the usual 

renormalizat ion-point mass \x of perturbation theory, which arises precisely because perturbation 

theory requires a cutoff. We believe that this parameter persists even in the homogeneous, finite 

Schwinger-Dyson equation because, as mentioned above, the categorical statement Z\ = is not 

true. There are other forms of the vertex Schwinger-Dyson equation in which it is essential to 

introduce a cutoff and cancel the cutoff dependence of Zi against that of various integrals involving 

the vertex. In other words, the homogeneous finite Schwinger-Dyson equation we use recognizes 

its heritage as a renormalizable, not super-renormalizable, theory. It is, of course, possible that 
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combining the study of the Schwinger-Dyson vertex with some other Schwinger- Dyson equation 
of the theory would fix the free parameter we find, but we have no evidence for that. 

A final result worth commenting on is the partial uncovering of instanton phenomena in sums 
of graphs — not just bare graphs, but dressed graphs, which leads to the incorporation of what 
would be called one-loop RG effects at our level of investigation. It is not straightforward to deal 
with instantons in the usual way, improving semi-classical results with the perturbative RG. For 
one thing, these perturbative corrections to instantons are not under control any more than they 
are in other sectors of the theory; there are IR renormalons which must be tamed by masses as in 
Eq. |2|. Furthermore, there are no exact instantons in the massive theory and perhaps most crucial, 
standard instanton techniques are really only applicable when all external momenta vanish ||. 
Ultimately, instanton effects at all momentum scales will have to be dealt with by Schwinger- 
Dyson methods. What we leave completely open here is how unitarity or some other physical 
effect provides a definition for the divergent sums associated with instantons. One would not be 
surprised to find a modification like Eq. [| entering for instantons as well as IR renormalons, at 
least in a well-behaved theory like QCD. 

Which of our results might persist in the theory of real interest to us, QCD? Even forgetting 
about confinement, the major complication which sets QCD apart from 0g is gauge invariance. 
Using the pinch technique @||, however, it is possible to define gauge-invariant proper vertices and 
self-energies which are related by naive (ghost-free) Ward identities. For these special vertices and 
self-energies Z\ = Z2, and both these renormalization constants vanish like (In Auv) -1 ^ 2 as Auv — ► 
00; one can, in principle at least, write down finite homogeneous Schwinger-Dyson equations for the 
gauge-invariant vertex corresponding to the running coupling, and these should show phenomena 
similar to what we did for 0|. For example, generation of a QCD mass scale is expected to yield 
a renormalon cancellation mechanism like that expressed in Eq. ||] for Borel-transform renormalon 
poles; indeed, this mechanism has already been invoked [|K|] on other grounds. It is also reasonable 
to expect that we can find instanton phenomena in graphical sums, just as for 

As for persistence of a free parameter like fi, the situation is not so clear. The vertex and 
propagator are not independent as they are in 03, which might lead to further constraints. One 
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expects that QCD (without quarks) has only one free parameter, the RG mass Arg, i n terms of 
which all physical quantities are determined. This does not seem to be the case for </>g, which has 
both and the mass M free, at least at the level of our investigation. 



II. PRELIMINARIES 

In this section we make some general remarks about the nature of the perturbation series 
for g<p 3 theories. We encounter non-Borel-summable series in g of the usual asymptotically free 
type, but in certain cases these series define one or more entire functions of g. We also set the 
stage for the investigation of 0|, discussing some aspects of the Schwinger- Dyson equations for the 
proper vertex, the proper self-energy, and a vertex formed from these which is RG invariant and 
corresponds to the running charge. 



A. General 6 Theories 



Theories of 3 type are well-defined, in general, only for purely imaginary coupling where they 



sometimes make sense physically [0)0] • ft * s sometimes possible to continue the theory to real g 
by imposing definitions (typically, the order in which integrals are to be done) which may or may 
not make sense physically. Consider as an example a zero-dimensional model with two fields x 
and y with a partition function given by 

/•OO 9 9 9 

Z(g)= / dxdye~ x y+9xy (3) 



which is well-defined if g is imaginary. We can give two forms to Z depending on which order the 
integrals are done; integrating over y first gives 

Z(g) = V^ f°° dx e - 2+ s 2 - 4 /4 (4) 

J — oo 

which is a wrong-sign 4 theory. Expanding Eq. |] in powers of g 2 gives rise to a typical asymp- 
totically free series with terms behaving like N\ N b (ag 2 ) N where a and b are fixed constants with 
a > 0. It is not obvious how such a divergent series is to be defined. 

If we instead perform the integral over x first in Eq. ^ we get 
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Z(g) = V^ / dy(l - gy)- 1 / 2 e^ 2 . (5) 

This appears to be singular at y — 1/g (the classical saddle point) of Eq. |3] but in fact Z(g) as 
defined by Eq. |5] is an entire function of g: as g changes from a pure imaginary value to the real 
axis, or anywhere else in the complex g-plane, the contour can be deformed to avoid the potential 
singularity. Actually, Eq. |5] defines two entire functions, depending on whether one begins with g 
on the positive or negative imaginary axis. One or the other of these possibilities may be singled 
out by physical considerations, or a linear combination may be used (e.g., the average of the two 
entire functions is real for real g). 

These remarks can be generalized to d — 6 theories of the type g<p\ip\ 2 where is a single- 
component Hermitian field, and ip is a single-component complex field. The Euclidean partition 
function is 



J (d(f)dipdip) exp — J d e x ^| d^ip \ 2 + -(d^<j)) 2 + g<f)\ip\ 2 + mass terms 



(6) 



The functional integral over is free, and when done reveals a wrong-sign non-local theory 
for ip, analogous to Eq. [5|. But doing the (free-field) functional integral over ip gives 

„o 1 



J (d(f>) exp - J d 6 x (-(<9 M 0) 2 + -M 2 2 - Tr ln(-D + m 2 + #)) 



(7) 



Superficially this would appear to define an entire function of g, since the quadratic terms in 
the action dominate the logarithm. However, the argument of the logarithm may vanish (when 
<fi ~ M 2 /g) and no firm conclusion can be drawn. Moreover, dimensional transmutation tells us 
that g must disappear, to be replaced by a running coupling. Still, it is intriguing to speculate that 
asymptotically free theories in general may have matrix elements and Green's functions which are 
entire functions of the running charge. If some such speculation is correct, it can only be developed 
by methods completely divorced from perturbation theory. 



B. The Schwinger-Dyson Equations for 0g 



Throughout this paper (except for section IV) we work only at Euclidean momenta; the action 
for our theory is 
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s 



/^[i(^) 2 + ^-0 2 + |0' 
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(8) 



We initially focus our attention on the Schwinger- Dyson equations for the renormalized propagator 
A and proper vertex T. Later we will construct the running coupling from a particular combination 
of these quantities. The equation for T can be written in several equivalent ways. In Fig. 1 we 
express T in terms of a four-point function which is two-particle-irreducible in the k\ channel. The 
four-point function can then be decomposed in terms of skeleton graphs and three-point functions. 
In this and other figures, a small blob denotes T, a line denotes A (both of which are renormalized 
quantities) and Z\ is the vertex renormalization constant. Figure 2 shows the analogous equations 
for T in terms of the one-particle-irreducible four-point function and Fig. 3 shows the equation for 
the inverse propagator. 

One of our major concerns will be that Z\ vanishes in an asymptotically free theory. This is 
not entirely an elementary matter since uncritically setting Zx — in Fig. 2 or the first equation 
for A -1 in Fig. 3 appears to give zero. The resolution, of course, is that the theory is defined first 
with an ultraviolet cutoff Auv which renders Z\ finite. In reality, as Auv — > oo the product of Z\ 
with a quantity involving a single bare vertex actually remains finite since the Auv dependence of 
such quantities cancels that of Z\. It is permitted, however, to set Z\ = in those forms of the 
Schwinger-Dyson equations in which no bare vertex appears (e.g., for T in Fig. 1 and the second 
form for A -1 in Fig. 3). The reason is that in an asymptotically free theory the remaining dressed 
graphs of the Schwinger-Dyson equation are all individually finite. This does not mean that the 
Schwinger-Dyson equations are unambiguously finite, since the Schwinger-Dyson equation for T, 
expressed in terms of T and A, has infinitely many terms and the sum actually diverges, as we 
discuss in Section IV. 

Let us concentrate for the moment on the lowest-order Schwinger-Dyson equations for T and 
A -1 , shown in Fig. 4, where we have now set Z\ = 0. Introduce the notation 



where M 2 is the renormalized mass. We anticipate, guided by one-loop RG-improved perturbation 



A-^A; 2 ) = {k 2 + M 2 )Z 2 {k 2 ) 



(9) 



theory ||, that the solutions to the Schwinger-Dyson equations of Fig. 4 behave in the UV (k 2 3> 
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M 2 ) as various powers of a logarithmic function D(k 2 ): 

D(k 2 ) = 1 + bg 2 \n(k 2 /M 2 ) (10) 
T -> LT Q , Z 2 -> IT 7 . (11) 

In Eq. |1(], b is a positive number given below; it is the same coefficient as in the one-loop /3-function 
(3 = —bg 3 + ■ • -. Any finite mass may be used to set the argument of the logarithm in Eq. |TU| and 
in consequence the 1 on the right hand side of this equation is irrelevant; we will drop it. 

Insert these forms into the equation for T of Fig. 4a using a simplified form suitable for finding 
the UV behavior: 

r{k2) = n r^ ( ^ inp2)37 ' 3a = ^ 2infc2 )" a • ( i2 ) 

Doing the integral reveals consistency, provided that 

1- 2a + 37 = 0, (13) 
2(4n) 3 ab = 1 . (14) 

Note that if the condition of Eq. [13] is satisfied, g 2 drops out of the equation so we might as well 
drop it in D of Eq. |10|, along with the 1 previously dropped. That is, there really is no coupling 
constant in an asymptotically free theory, only mass scales. 

We need one more equation to determine a, 7, and b, which is furnished by the A -1 equation. 
This is slightly ticklish because it has a quadratic divergence, subsumed by mass renormalization. 
In a g4>\ip\ 2 theory one may avoid direct consideration of the A -1 equation by coupling ip to the 
electromagnetic field in the usual way, computing the lowest-order electromagnetic vertex shown 
in Fig. 5, using the Ward identity 

q^\q,p,p > ) = A- 1 (p)-A- 1 (p') (15) 

and then setting p' 2 = M 2 , so A _1 (p /2 ) = 0. It turns out (see below) that Z 2 = also, so the 

equation is homogeneous of Baker- Johnson- Willey type. We leave it to the reader to conduct an 

analysis like that which led to Eq. [12|; it gives the same result as Eq. [16| below. 

Either as sketched above or by analysis of the A -1 equation in Fig. 3 one finds a new condition 
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(47r) 3 7 fo = 1/12 . (16) 

Then one finds 

a = 2/3, 7 = 1/9, b = ^ (17) 

which are the usual results || of one-loop RG-improved perturbation theory, with /3-function 
— bg 3 + • • •. These constants also appear in the cutoff dependence of the perturbative renormal- 
ization constants 

Z 1= (l + bg 2 H^/^)y 2/ \ Z 2 =(l + bgHn(Al Y /^))' 1/9 (18) 

both of which vanish as the cutoff Auv ~~ > oo. 

At this point it is convenient to introduce a special combination of T and Z 2 , corresponding 
to the running charge 

-a / u\- T(k 1 ,k 2 ,k 3 ) , . 

g{k 1 ,k 2 ,k 3 ) = g (19) 

y /Z 2 (k 1 )Z 2 (k 2 )Z 2 (k 3 ) 

We keep the explicit g only to observe that, with D = bg 2 In k 2 , T ~ D~ 2 ^ 3 , Z 2 ~ D^ 1 ^ 9 , all the g's 
cancel in Eq. [19]. We now find the asymptotic behavior, when all ki scale like a large momentum 

k, 

' 2 ~* MnJ? ' (20) 

This is the usual one-loop running coupling. If constructed in perturbation theory, g would be 
explicitly independent of the renormalization point \x. 

With Z\ — 0, the Schwinger-Dyson equation for g (the solid blob in Fig. 6) looks like that 
for T with the important difference that now the propagator lines correspond to free propagators 
with the renormalized mass. This seems to lead to a quite remarkable circumstance in which the 
Schwinger-Dyson equation for g is completely self-contained. However, this is misleading since 
the infinite series of terms in Fig. 6 does not have a sum which is well-defined. In fact, this series 
contains the basic instanton phenomenon of 0g, as we will see in Section IV, but we will never 
need to invoke instantons explicitly. We will also see in Section IV that Schwinger-Dyson graphs 
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for g having 2N + 1 insertions of g behave at large momentum like (In k 2 )~ N ^ 2 , so that the leading 
UV behavior is correctly captured in the cubic graph. 

In the next section we will investigate toy models of the Schwinger-Dyson equation somewhat 
similar to the integral of Eq. [12]. If we uncritically write such an equation for g with free propagators 
(7 = 0) and test whether g = (blnk 2 )^ 1 ^ 2 is a solution, we find that everything works except that 
the integral is 4/3 larger than it should be. This is because the simple form assumed for g does not 
distinguish the momenta on the legs of g which are in general different; this distinction matters 
for g with its separate powers of Z 2 {ki)~ l l 2 , not all of which involve momenta which are integrated 
over. In the toy models of the next section we will compensate for this by multiplying the right 
hand side of integrals cubic in g by 3/4. 

Why do we need toy models at all? The answer is that even without the spin complications 
of QCD, and saving only cubic nonlinearities, the full Schwinger-Dyson equation for (f>\ is quite 
complicated. Here we set up this cubic equation schematically, before going on to the toy models 
in the next section. 

The proper vertex T(ki, k 2 , k 3 ) depends on three scalar variables k\,k\,k 2 in a way which is 
constrained by causality and spectrum conditions. These may be enforced with the Nakanishi 



representation |13| , which is derived from the Feynman-parameter representation of graphs but 
which is expected to hold even for non-perturbative processes. For the vertex it reads: 

r(h, k 2 , k 3 ) = [ dz 1 dz 2 dz 3 dcr5(l - Z\- z 2 - z 3 ) — , —ri~r~TTT — H ■ ( 21 ) 
J a + Z\k{ + z 2 k 2 + z 3 k 3 

The spectrum is expressed through the support of h. In a massive theory, the lower limit of the 

a integral depends on masses and Zi whereas the upper limit is infinity. The asymptotic behavior 

T ~ (ln/c 2 ) -2 / 3 requires that h ~ (In cr) -5 / 3 at large a (in perturbation theory h — > constant, and 

Eq. needs the usual subtraction). One then uses Eq. and the Lehmann representation for 

the propagator to write out the Schwinger-Dyson equation of Fig. 6. If the loop momentum is 

called p, the momentum integral depends on three scalars (p 2 ,p- k\,p- k 2 ). One might then try to 

solve the equations numerically, but it is a substantial project. If each k 2 takes on 100 values, one 

has to do 10 6 three-dimensional integrals numerically for every iteration. We have not attempted 

this, but have studied some toy models which we are convinced have the qualitative features as 

11 



the real (cubic) Schwinger-Dyson equation. We now turn to those models. 



III. THE TOY MODELS 

We have examined several models based on simplifying the cubic Schwinger-Dyson equation of 
Fig. 1; we will describe a few here. A simplification common to all the models is to assume that 
r(&4, &2, ^3) depends only on k\ + k\ + k\ which corresponds to fixing all Z{ = 1/3 in the Nakanishi 
representation of Eq. ^T. 

In this section we restrict ourselves to the special momentum configuration k\ = 0, k 2 = — k^ = 
k, and we will be interested in the whole (Euclidean) range < k 2 < 00. We will always use 
free propagators in the model Schwinger-Dyson equation, which can either be interpreted as an 
approximation to the propagators in the Schwinger-Dyson equation for the proper vertex T, or 
as a special kind of simplification of the Schwinger-Dyson equation for the running coupling g, 
where free propagators are required (see the discussion around Eq. |T9|) , but where the assumption 
of dependence of g only on ^k 2 cannot really be correct because of the (Z2(/ci)Z 2 (/c 2 )Z 2 (/c3)) -1 / 2 
factor in its definition in Eq. [19|. We choose to present our results as referring to g; the error 
committed in assuming g depends only on kj is quantitative but does not affect our major 
qualitative conclusions. The reason for our choice of interpretation is that all model solutions 
behave like g in the UV region, that is, like (lnfc 2 ) -1 / 2 . 

Various models differ from each other by further simplification of the momenta appearing in 
internal vertices. The Schwinger-Dyson equation with correct momentum assignments is shown 
in Fig. 7, and corresponds to the integral equation 



(Pp G(2(p + kf)G 2 (p 2 + k 2 + (p + k ) 2 ^ 
7r3" (Go+fc) 2 + m 2 ) 2 (p 2 + M 2 ) 



G{ 2k 2 ) = I ^ 1 ;; 9 y Z\ J„ ■ ( 22 ) 



where the normalization is chosen so that G — > (In A; 2 ) l l 2 for k 2 3> M 2 ,m 2 so that we interpret 
G as b +1 ^ 2 g. Later we will present numerical solutions of Eq. |22, but for now it is much more 



important to have an analytic understanding of simplified versions of this equation. The reason 
for using two masses M and m will be made clear later. 
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A. The Massless Model. 



The first model to be considered is the massless model (M — m — 0) where, to make things 
more tractable, we set k = in the arguments of G on the right hand side of Eq. [22| to get 

It can be shown, by an uninteresting argument bounding G from above and below, that the 
solutions to Eq. |22] and Eq. ^ have the same behavior for large k. 

Performing the angular integrals in Eq. |23| and introducing the variables 

t = ln(fc 2 //i 2 ) t' = Mp 2 /n 2 ) (24) 

where \i is an arbitrary mass scale, we obtain the one-dimensional integral equation 

1 r°° 1 rt , 

G(t) = - dt'G 3 (t') + - dt'e^-^G^t'). (25) 

2 J — oo 2 J — oo 

Differentiating twice, one finds that G obeys the equation 

G + 2G = -G 3 (26) 

where G = dG/ dt. This same equation could be obtained by noting in Eq. |23| that [— Air 3 (p+ A;) 4 ]" 1 
is the inverse of the d'Alembertian operator in six dimensions. 

We must point out that Eq. ^6] actually makes no sense if we require G(t) > (as physical 



considerations suggest). Perhaps the most intuitive way to see this is to think of Eq. as the 
differential equation governing the 'position' G of an anharmonic oscillator as a function of 'time' t: 
except for the trivial solution (G = 0) all solutions to this equation inevitably develop unbounded 
oscillations as t — > — oo. The problematic region t — > — oo corresponds to an IR divergence at k = 0. 
It is this divergence which is responsible for the renormalon behavior we will shortly encounter. 
The t — > +oo solutions are candidates for the large-fc behavior for the massive equation of Eq. [22|. 

There are at least two fundamental types of large-t solutions to Eq. and another which is 
a hybrid of the first two. In the first type, the G term in Eq. ^6] is treated as a perturbation to 
the equation 2G = —G 3 which has a solution G = t^ 1 ^ 2 = (ln(/c 2 / /j 2 ))^ 1 ^ 2 , as anticipated from 
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field-theoretic perturbation theory. Including the G term as a perturbation, we exhibit several 
more terms of this solution, 

G = D- 1 / 2 ] D = t + ^ Indt-^l- ^lnpat)+... . (27) 

Here the q are constants not determined by consideration of large t alone. 

The behavior of Eq. |27| is, except for the numerical coefficients, precisely what the RG and 
perturbation theory would give for the running coupling. In fact, it is possible to find the (3- 
function for the vertex g by noting that 

(3(g) =2g = Ib^G (28) 



which leads immediately from Eq. 26 to a first-order differential equation for (3(g): 



(3(l + lf)=-bg 3 . (29) 



4dg, 

The power series solution beginning with bg 3 shows a non-Borel-summable factorial divergence 
associated with renormalons, 

/%)=-£w QV +1 . (30) 

where, for large N, = means equality up to fixed powers of N and overall multiplicative constants. 

We will return to this solution and the associated g shortly, showing how to set up Eq. ^ directly 

as an integral equation for a Borel transform. 

There is another class of solutions to Eq. |26| which shows no relation to perturbation theory; 

this class arises from by treating the — G 3 term in Eq. |26| as a perturbation. The solution reads 

A p- 2t A 2 A A A 6 p- 12t 

G = : Di = 1 + — e~ 4 * + e~ 8t + + . . . . 31 

D 1 ' 24 5760 967680 V ; 

where A is an arbitrary positive constant. Note that this solution behaves as inverse powers of k 4 , 
powers which one would like to associate with condensates and higher-twist terms. 

Finally, there is a hybrid solution combining the previous types of solutions of Eq. ^ and 
Eq.0, 

G = L>2~ 1/2 ; D 2 = t + - lncit + . . . + 2At 3 e~ 2t + ... , (32) 
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where again c\ and A are arbitrary. The solutions of the real Schwinger- Dyson equation also ought 
to have a combination of perturbative powers of logarithms and powers of k 2 . The existence of 
such hybrid solutions show that more numerical analysis is not very revealing, since it is hard to 
separate powers and exponentials of t (i.e., of ln/c 2 ). 

An interesting way of looking at our models is to introduce the Laplace transform of G(t): 

roo 

G(t) = / da F(a) e~ at . (33) 
Jo 

Since in the lowest order of approximation t ~ (bg 2 )~ l , the Laplace transform F(a) can also 
be interpreted as the Borel transform. At this order this appears to be a trivial matter, since 
G = b l l 2 g, t = (bg 2 )~ l leads to F(a) = (na)^ 1 ^ 2 . But as one goes on more interesting phenomena 
appear. 

By recognizing e~ at = (k 2 /fi 2 )~ a and using standard integration formulas, one can express the 
massless integral equation Eq. ^3] in the form 

01 r 

F(a) = « / \. d v\ F(ay 1 )F(ay 2 )F(ay- i ) (34) 

z — a J 

where [dy] = f x dyidy 2 dy 3 5(1 — J2Vi)- Here are Feynman parameters used to express the 
Laplace transform variables (aj = (xyi) of the three vertices on the right hand side of Eq. |23| . 
Perturbation theory corresponds to a < 1, in which instance we replace 2 — a by 2, and verify 
that F = [jia)^ 1 / 2 (i.e., G = t^ 1 / 2 ) is the small-a solution. But there clearly is some sort of 
singularity at a = 2, which actually causes Eq. ^ to be meaningless (just as integrals leading 
to IR renormalons are meaningless). Consider the "approximation" of replacing the F(ayi) by 
their perturbative values on the right hand side of Eq. the result is an F(a) behaving like 
a _1 / 2 /(2 — a). We can then construct the Laplace transform for (3(g) ~ G, and in e~ at replace t 
by (bg 2 ) -1 . The result is a Borel-transform representation of /3(g), 

roo q,1/2 

(3(g) ~ / da e-^ 2 . (35) 

Jo 2 — a 

The pole at a = 2 is, of course, what gives rise to the factorially-divergent behavior of (3 shown in 
Eq. 

A major theme of our work is the demonstration of how such Borel-transform singularities are 

cured in the massive theory, as we know they must be. So now we turn to massive models. 
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B. Massive Models. 



We now return to the massive equation Eq. If m 7^ it is no longer possible to find a local 
differential equation as we did for the massless case, and analytic progress beyond the leading UV 
behavior seems impossible. However, when m = we can convert Eq. |22| to a differential equation 
(or a corresponding integral equation in one variable) just as we did when both M and m were 
set to zero. As long as M 7^ 0, there is no IR singularity in Eq. |22], as one easily checks by looking 
at k — 0. We will, using the Laplace transform technique, translate this lack of IR divergence 
into a cancellation mechanism for the poles of the Laplace-transform kernel which appeared in the 
massless case (see Eq. |4j). 

We have actually looked at several variants of Eq. |22|, which differ in the choice of m and of 
the arguments of the G's under the integral. These are conveniently summarized via: 

g(2*Wj£ .^L. (36) 



tt s ((p + k) 2 + m 2 y (p 2 + M 2 ) 
The models we will discuss are, for M^O: 

A) m = 0, F = G 3 {2p 2 ) . 

B) m = M, F = G 3 (2p 2 ) . 

C) m = M, F = G (2(p + A;) 2 ) G 2 {p 2 + k 2 + {p + k) 2 )). 

The k— independence of F in models A and B gives rise to one-dimensional integral equations which 

we analyze both analytically and numerically. Model C, on the other hand, is Eq. [22| with m = M 

which gives rise to a two-dimensional integral equation which we attack numerically. In addition to 

models A-C we have also examined intermediate cases where, for example, F = G 3 (2(p + k) 2 ) with 

m = M or F = G 3 (2(p 2 + k 2 )) with m = M. These additional cases introduce a dependence 

into F which is still mild enough to yield one-dimensional integral equations. However, since 

these additional models have all the same qualitative features as models A-C but yield little extra 

insight, we will not discuss them further. 

Only model A will be discussed in much detail. One can show that all models lead to the same 

leading UV behavior (G — > (In A; 2 ) -1 / 2 ), and numerical work displayed later shows that the models 
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differ quantitatively but not qualitatively in the IR. All have the curious feature, mentioned in 
Section I, that the solutions to each model depend on a single real parameter. 

Consider now model A. Similar to our treatment of the massless model, we introduce the 
variables 

t = ln hvH' *= ln T (37 > 



where we have fixed the analog of the arbitrary mass scale fi in Eq. |24] to be equal to M. In terms 



of these variables, Eq. |36] for model A reads 

i r°° , , „ l /•* ,r/e 4 ' - in 2 



G(t) = -J dt'G 3 (t') + -J o dt> 



2 Jo 2 Jo 

which is equivalent to the differential equation 



e* - 1 



G 3 {f) . (38) 



(1 - e~*)G + (2 + e- l )G = -G 3 (39) 

with the boundary conditions 

1 f°° 

G(0) = - / dt'G 3 (t'), G(oo) = . (40) 

2 jo 

This differential equation differs from the massless equation Eq. [26| by the appearance of e~* in 
the coefficients; it can, of course, be derived directly from Eq. |38], just as we did for the massless 
equation Eq. 

It is not hard to see that G(oo) = is not really a boundary condition since all solutions of 
Eq. [39| vanish at t — oo. Moreover, Eq. ^0] does not constrain G(0); as long as G(0) > 0, there is a 
solution to Eq. |39] satisfying Eq. ^0|. So G(0) itself can be chosen as the parameter of the solution. 
Once G(0) is picked, the solution is unique. All of the models listed above have this feature. 

To see how this parameter arises for Eq. [38|, consider a slightly different equation with an 
explicit parameter a: 

G(t, a) = a-- f dt'G(t', a) 3 + - f dt' ( ^— ^ ^(f, a) 3 . (41) 



2 Jo ~" K ' ' '2 Jo Ve*-1 

For t 1 the factor (e*— l) 2 / (e*— l) 2 is exponentially peaked ait' = t so it is a good approximation 

to evaluate the relatively slowly varying G(t', a) 3 at t' = t and pull it outside the integral to get 
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1 /■*>! 1 ~ 

G(t> l,a) = a-- / rft'G(t',a) 3 + -G(t > l,a) 3 . (42) 

2 jo 8 



If we use the fact that G(t 3> 1, a) ~ and let t go to infinity, we find 



a = - / dt'G(t', a) 3 . (43) 
2 jo 

In other words, if G(t, a) is a solution to Eq. [y] then G(t, a) automatically satisfies Eq. [38]. 

We do not know how to solve Eq. [39] analytically, so yet another presentation of its contents 
is needed to understand the taming of IR renormalons. We proceed as in the massless case 
by introducing the Laplace transform Eq. [J3], with t given by Eq. |37|. Introducing a Feynman 
parameter x for the p 2 + M 2 propagator and performing the momentum-space integral, one gets 



G(t) = da a [dy]F(ay 1 )F(ay 2 )F(ay 3 ) / dx — — — — . (44) 

jo J Jo (1 + x{e l — I)) 

The x-integral is elementary, and gives: 



(45) 

This is essentially the massive version of Eq. for F(a), which is singular at a = 2. But 
there is no singularity in the integrand of Eq. [E| for any positive a. We now see the qualitative 
(even though quantitatively inaccurate) solution to the IR renormalon problem, expressed as the 
cancellation by what amounts to condensate terms of poles in the Laplace transform. It is evident 
that the same sort of cancellation takes place in the /3-function, replacing e~ a ^ bg2 — e~ 2 ^ bg2 (and, 
given Eq. |45|, another term is added involving a cancelled pole at a = 1). 

Note, by the way, that it is elementary to derive the differential equation Eq. [39] from Eq. |45|: 
Just form 

One can look for solutions in the IR or UV by the same method as for the massless equation. There 

are always well-behaved power-series solutions in t, with G(0) as an undetermined parameter; 

these are useful in the IR. In the UV, there is the same solution in powers of t~ l l 2 (modulo 

18 



logarithms) as for the massless case Eq. ^7], but with exponential corrections. There is also an 
exponentially-vanishing solution based on treating the — G 3 term as a perturbation; the leading 
term is (e* — 1)~ 2 (cf. the massless case Eq. |31]). But it appears from numerical studies that this 
exponentially vanishing solution develops a singularity for some t > 0, which we exclude. Thus 
the only possibility is the usual UV solution of Eq. |2?|, at least for generic values of the vertex 
momenta fcj. 

It might happen that <fi\ has a bound state with the quantum numbers of the 0-field, for some 
particular value of k\ ^ 0. In that case, the Schwinger- Dyson equation for the bound-state wave 
function might decrease exponentially in the UV; we have not studied this possibility in detail. 

C. Numerical Solutions. 

We have used numerical techniques to solve the integral equations associated with the three 
models A-C discussed above. Of these only model C gives rise to a two-dimensional integral 
equation; the others can be put in the one- dimensional form 

G(t) = / dt'K(t,t')G 3 (t') (47) 
Jo 

with different kernels K. To solve Eq. f|7|, we put it in the form 

roc _ 

Git) = a+ dt' K(t,t')G 3 (t') (48) 
Jo 

where K = K — 1/2 and we treat a as a free parameter, as in Eq. |41]. For all kernels, K falls 
off rapidly when t ^> t' and one shows a = | J °° dt G 3 , as before. The two dimensional integral 
equation associated with model C is treated in a completely analogous manner, with a similar 
definition of a. We have explored thoroughly the range < a < 2, and we will show results for 
a = 1,2. All results were calculated to 16 digits precision. For a > 2 it takes progressively longer 
for our numerical routines to converge so we cannot exclude the possibility of an upper bound to 
a above which solutions cease to exist. 

Fig. 8 compares models A (m = 0) and B (m = M). For model B G(0) = a (see Eq. ^0[) which 
is larger than G(0) for model A, increasingly so as a gets larger. This comparison shows the effects 
of varying mass assignments. 
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Fig. 9 compares models B and C. which have all masses set to the same non-zero value but 
differ in the arguments of the G"s in the Schwinger-Dyson equation. For the same a value (1 in 
this case) the differences are remarkably small throughout the entire range of momenta. Results 
such as these give us confidence that toy models, such as model A and B which can be studied 
analytically to some extent are likely to be representative of the full Schwinger-Dyson equation, 
even at zero momenta. 

The next question to ask is what happens when one goes beyond the cubic term in the 
Schwinger-Dyson equation. 

IV. AN APPROACH TO SCHWINGER-DYSON VERTEX GRAPHS OF ALL ORDERS 

In this section we take the first incomplete steps towards an analysis of graphs of all orders in 
the Schwinger-Dyson equation for g or for the proper vertex T. These Schwinger-Dyson graphs are 
expressed by writing all perturbative vertex graphs which have no vertex or self-energy insertions; 
a few such graphs are shown in Fig. 6. Then in the full Schwinger-Dyson equation each such 
graph is replaced by one of the same topology but with full vertices and propagators. In what 
follows we do not need to consider the cubic graph which was the concern of previous sections. 

All remaining vertex skeleton graphs have the property that the minimum number of lines in 
any closed loop is at least four; graph theorists say that such graphs have girth four. It turns 
out that such bare skeleton graphs, in perturbative 0g, have exactly one UV logarithm (lnAyy, 
where Ayv is a UV cutoff) and no IR logarithm (InM 2 ). This follows straightforwardly from an 
elementary remark about (p\ graphs, as we now show. 

Every perturbative skeleton graph of N loops, (3N lines, 2N + 1 vertices) has the following 
Feynman-parameter form in Minkowski space, which is for the moment convenient: 



where [dx] = 5(1 — Y. x i) I\i=i d x i- U is the graph's determinant (a positive sum of monomials of 
order N in the Xj), and <fi (not to be confused with the field 0) is a sum of positive polynomials in 
the Xi times scalar momentum functions k a ■ kf,. The single UV divergence is explicit, and can be 




(49) 
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removed by taking the imaginary part of T° N , replacing the logarithm by a ^-function. Any other 
divergences in the imaginary part must come from the vanishing of U. But, as is well known, U 
can vanish only when at least all the parameters of a loop in the graph go to zero. Let {xj} be the 
parameters of one and only one such loop of n lines, and replace Xj by Xxi, X —>■ 0. Then U vanishes 
like A. But [dx] behaves like X n ~ 1 d\ times irrelevant factors, and the integral is finite if n > 4, 
that is, for girth-four graphs. If all the parameters of two complete (possibly overlapping loops) 
are scaled with parameters Ai, A2, then U vanishes like A1A2, and the argument can be repeated, 
with some necessary variations. The reader may easily check that only for skeleton graphs the 
total number of distinct lines in i complete loops is at least 3^+1 for any skeleton graph of order 
N > i. Then the one-loop scaling argument given above can be repeated, and one shows there 
are no singularities in the Feynman parameter integral. (The case N = t is trivial, since not all 
the Xi can be scaled to zero.) 

The Minkowski-space imaginary part of T° N in Eq. 4£ is, when all masses are equal to M, 

lmT% = constant x g m J ^ - M 2 ) . (50) 

Consider now the kinematics of Fig. 1 where k\, — k 2 , and — k$ are forward timelike, with k\ = 
s > AM 2 , k\ — k\ — M 2 . One easily shows k\ ■ k 2 = kx ■ k 3 — s/2, k 2 ■ k 3 = s/2 — M 2 , which is 
enough to show that the argument of the ^-function is positive (trivially so at M = 0) and thus 
the ^-function can be replaced by unity. We then reduce the problem of finding Imr^ to summing 
the integrals 

[dx] 



v? < 51 > 

where i runs over all the N — th order skeleton graphs. This problem (with U^ 3 — > U^ 2 ) has been 
addressed for <p\ [0, using the following steps which we will repeat here. 

1. Count all the N th order graphs, and note that the number of these which are not skeleton 
graphs vanishes at large iV relative to the number of skeleton graphs. We denote this number as 
Qn- 

2. Find, by a combination of analysis and numerical work, the number C, of monomials in Ui. 

This number is called the complexity of the graph, and is equal to the number of spanning trees. 
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3. Show that 



_ J[dx]U r _ (3iV-l)! 
{U)t = JJdxT ~ C \AN-1)\ ■ (52) 

That is, every monomial in Ui contributes exactly the same to (U){. 

4. Use a Holder inequality to show 

hf.-jU)!' < 53 > 

5. Observe that the distribution of G\ for all the iV-loop graphs centers around an average 
value (C). 

We then have, inserting the usual powers of 2it, 

lmT ^(wvT QN wr (54) 



Now Bender and Canfield JTj] have shown, in a deep graph-theoretic result, that 

/On N 

QN=(n) N\ (55) 



,2 

and we have made a numerical analysis of large- N graphs to find 

/16\ N 

(C)=( T ) . (56) 

(Here A N = B N if A N , B N are each of the form c D (A^!) Cl c 2 JV A^ C3 (l + 0{l/N)) at large N, and the 
constants c\ and c 2 are the same for and B^.) 
Putting everything together, 

N 2 11 



Imr^ > A^! 



a g 



— ~ 0.312 . (57) 



L(4tt) 3 . 

As expected for an asymptotically free theory, this is not Borel-summable. 

Even though extracted from a high-energy process in Minkowski space, the result of Eq. [57] 
looks exactly like what one would expect from a lowest-order (semi-classical) instanton calculation 
which requires zero external momenta and masses. There is no sign of external momenta or masses 
in Eq. [57], which renders this formula compatible with the Lipatov procedure. So let us compare 
the coefficient a in Eq. [57] to the instanton result. As is well known 0, there are instantons in 



massless 0g, of the form 



22 



48p 2 



(58) 



g[(x - a) 2 + p 2 } 2 



whose action is 



\5j g 2 



(59) 



This means that using the Lipatov procedure one would find an expression of the form of Eq. [57] 
for Imr^, with a = 5/12 = 0.417. Our lower-limit value 2 11 /3 8 is about 3/4 of this (presumably 
correct) value, and further numerical work, which we have not done, would bring this lower-limit 
value closer to the true value of a. 

The fact that there are no logarithms in Imr^ (or only one in T^) is of vital importance for 
analyzing the Schwinger- Dyson equation, when full propagators and vertices are used to dress 
the bare skeleton graph. This amounts to inserting 2N + 1 factors of g (see Eq. [19]), where 
every g behaves like (b In k 2 )~ x l 2 for the appropriate momentum k. One should then not be 
surprised to find that the iV-loop full skeleton graphs behave, for large external momentum k, like 
(lnfc 2 )-2( 2JV+1 ) +1 and that the sum of full skeleton graphs is of the form Eq. [57], with g 2 replaced 
by the (one- loop) running charge g 1 . The Schwinger-Dyson equation thus yields the results of 
RG-improved instanton calculations, with no explicit introduction of instantons. 

We will not show the derivation of this result in detail, but merely sketch it. As in Section III 
we make the approximation that every g depends only on k 2 , where ki are the vertex momenta, 
and introduce the Laplace transform 



Anticipating that the leading UV behavior of the N th order graphs is dominated by the UV 
behavior of the cubic graph studied in Section 3, we choose 



As before, the factor exp(— a In J2 kf) is written as (Y^kj) a and the momentum integrals are 
done after Feynman parameterization. A tedious analysis, not given here, shows that one need not 




(60) 



F{a) = (Trba) 



1/2 



(61) 



corresponding to g — > (b In k 2 ) l l 2 . 
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introduce 2N + 1 new Feynman parameters for the new denominators J^kf] instead, one shows 
that the correct asymptotic behavior is also found by replacing such denominators by those of 
neighboring propagators. The result is that many of the propagators in a bare skeleton graph 
appear to powers — (1 + cKj), where a, is the Laplace transform variable of a neighboring vertex. 
Finally one finds, setting the mass M to zero, and going to asymptotically-large k as in Section 3, 



9n = Q 



N 



- 1 ■ 


7 


.7r(47r) 3 6_ 





Ha~ ll2 d ai Y{Y, «i) [nr(l + Oi)]' 



t/ 3 



U 2 17 



(62) 



where <fi and C/ are the appropriate functions for the bare skeleton graph, and an irrelevant constant 
factor has been dropped. In the result of Eq. |UT|, Qn is the number of iV th -order skeleton graphs, 
given in Eq. 

In the kinematics of Fig. 7, is of the form of k 2 times a function of t/i. We scale the aij by 
on = ay u with Y,Vi = 1, and get 



where 



i. 

2' 



9N = J 


poo 

' da — 




o r 


- 1 ■ 


N . 


.7r(47r) 3 6_ 


Qn J 



a 



r(iv-i) 



e- alnk2 H(a) 



(63) 



Q N J[dy]Uy; 1/2 T(l + a) J [ -^Ux 



\k 2 U 



(64) 



In Eq. ^3| and Eq. Q we wrote T(J2 a«) = T(a) as a _1 r(l + a) and removed a factor T(iV — |) for 
simplicity. This allows us to assert that H(a) is regular at a = 0, and in fact one checks easily 
that 



ff(0) 



A' 







jV 



T7 3 " 



.7r(47r) 3 6 

which is just the function appearing in the sum of bare skeleton graphs. It is clear from Eq. 
that the leading asymptotic behavior of comes from a ~ 0, which gives 



(65) 



^jv = ^(0) J da- 



3 

a" 2 



-a In k 



= N\ 



(47r) 3 61nP 



(66) 



T(iV-|) 

using the previous result for H(0). 

This is as far as we will carry the analysis. We have no particularly good ideas for extending 

the toy models of Section III to all orders in a controllable way. Nor will we attempt to find some 
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physical way to sum the series of Eq. ^6] (based, e.g., on decay of the meta-stable vacuum). If the 
sum of terms is represented by a (formal) Borel transform, there arises an ambiguity of the form 
(k 2 )~ c , where c = (47i) 3 b/a. This is of the usual type associated with RG-improved instantons. 

V. CONCLUSIONS 

In the long run, perturbative QCD must fail because of its IR singularities, and some sort 
of new perturbative analysis of Schwinger-Dyson equations will have to be seriously attempted. 
This is much more difficult than 0g, because of the complications of spin, gauge dependence, and 
confinement. Still, we believe that some of our results will hold for QCD, in particular, the cancel- 
lation of poles in Borel transforms as in Eq. |2|. Currently there is considerable phenomenological 
interest || in renormalon problems, and it would be valuable to test the mechanism we identify 
as relevant, versus other prescriptions such as principal part integration. 

A similar mechanism may well hold for Borel transforms associated with instantons, although 
we are unable to say in 0g because of its instability. In QCD (or a similarly-stable asymptoti- 
cally free theory, like the d = 2 nonlinear a model), it is quite important to learn how to deal 
with phenomena canonically associated with instantons (e.g., the Lipatov technique) by other 
means than the usual RG-corrected semi-classical results, because the perturbative corrections to 
semi-classical instanton physics involve IR singularities such as renormalons. Understanding such 
phenomena directly from the Schwinger-Dyson equations is a worthwhile thing to do. Moreover, at 
the moment Schwinger-Dyson equations are the most straightforward way of dealing |7|] with the 
non-Borel-summable divergences at large Minkowski-space external momenta. We hope to report 
in the future on future progress in this direction, going well beyond the first steps of Section IV. 

We have seen that all the toy models of the Schwinger-Dyson equation have solutions with a 
free parameter, equivalent to the renormalization mass fi. Since fi comes in as a tool in dealing 
with amplitudes needing regularization, and since asymptotically free Schwinger-Dyson equations 
in their canonical form need no such regularization, this is a bit of a surprise. On the other hand, 
one would be surprised if a single Schwinger-Dyson equation — in our case, that for the running 
charge defined in Eq. [TP] — had a single solution for the three-point function, making no reference 
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to higher-point functions as the usual hierarchy argument specifies. It remains to be seen how this 
is settled by consideration of further Schwinger-Dyson equations. 
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FIGURE CAPTIONS 



1. Bethe-Salpeter equation for the renormalized three-point function (small blobs) ex- 
pressed in terms of the two-particle-irreducible (in the k\ channel) four-point function 
(large blob) and the vertex renormalization constant Z\. All propagators are fully 
dressed. 

2. Bethe-Salpeter equation for the renormalized three-point function (small blobs) ex- 
pressed in terms of the one-particle-irreducible (in the k\ channel) four-point function 
(large blob) and the vertex renormalization constant Z\. All propagators are fully 
dressed. 

3. Bethe-Salpeter equation for the renormalized inverse propagator in terms of the renor- 
malized three-point-function and the vertex renormalization constant Z\. All propa- 
gators are fully dressed. 

4. Lowest order approximations to the Schwinger-Dyson equations for a) the three point 
function and b) the inverse propagator. 

5. Cubic Schwinger-Dyson equation in g\ip\ 2 <p theory. 

6. Schwinger-Dyson equation for the running charge (solid blob). All propagator lines 
are free propagators with renormalized masses. 

7. Momentum assignments for the truncated Schwinger-Dyson equation for the running 
charge. 

8. Numerical solutions of the toy- model Schwinger-Dyson equation for models A and B 
for a = 1, 2 where a = 1/2 J dt'G 3 (f). 

9. Numerical solutions of the toy- model Schwinger-Dyson equation for models B and C 
for a = 1. 
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